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Abstract: We present a top-down string theory holographic model of strongly interacting 
relativistic 2 + l-dimensional fermions, paying careful attention to the discrete symmetries 
of parity and time reversal invariance. Our construction is based on probe DT-branes 
in AdS^ X S^, stabilized by internal fluxes. We find three solutions, a parity and time 
reversal invariant conformal field theory which can be viewed as a particular deformation 
of Coulomb interacting graphene, a parity and time reversal violating but gapless field 
theory and a system with a parity and time reversal violating charge gap. We show that 
the Chern-Simons-like electric response function, which is generated perturbatively at one- 
loop order by parity violating fermions and which is protected by a no-renormalization 
theorem at orders beyond one loop, indeed appears with the correctly quantized coefficient 
in the charge gapped theory. In the gapless parity violating solution, the Chern-Simons 
response function obtains quantum corrections which we compute in the holographic theory. 
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1 Introduction and Summary 

The AdS/CFT correspondence [1] has recently been used extensively as a tool to exam- 
ine the behavior of quantum systems at strong coupling [2]. Emergent relativistic 2+1- 
dimensional fermions in condensed matter systems, such as graphene [3] [4], topological 
insulators [5] , the D-wave state of high Tc superconductors [6] and optical lattices [7] lead 
to interesting dynamical questions, often in the strong coupling regime. 

For example, a model of the low energy (< lev) physics of graphene has massless 
relativistic fermions interacting with a Coulomb potential and emergent SU{A) symmetry, 

4 
S= j d'xY.i^k [l\idt - At) + VFi ■ (iV - ^)] V'fc - ^ y" d'x F.^-^F"^ (1.1) 

where we have set h = 1 = c. Electron dynamics are confined to the 2-dimensional sheet 
of carbon atoms which composes graphene. The nonlocal gauge field action arises from 
the fact that the photon propagates in the 3-dimensional space surrounding the graphene 
sheet, This model does not have SO(2, 1) Lorentz symmetry since the speeds of the massless 
electron and the photon are not identical, in fact, vp ~ c/300. To a first approximation, 
the Coulomb interaction is instantaneous and the magnetic interactions are suppressed by 
factors of ^. 

The field theory with action (1.1) is scale invariant at the classical level and it is 
known that, at least up to two loop order, it is a renormalizable quantum field theory 
in the sense that no counterterms with new interactions have to be introduced to cancel 
ultraviolet divergences [8] . Upon renormalization, the graphene speed of light vp requires a 
logarithmically divergent counterterm and it therefore becomes a scale-dependent running 



coupling constant which increases to larger values in the infrared limit. This running 
should be cutoff when vp reaches the speed of light. The resulting Lorentz invariant theory 
is conjectured to be a conformal field theory for any value of the dimensionless coupling 
constant, e. This is confirmed to two loop order in perturbation theory. 

However, to describe realistic physics of graphene, (1.1) is nominally a strongly inter- 
acting model. The strength of Coulomb interactions are governed by the graphene fine 
structure constant, which (if we use the vacuum dielectric constant eo = 1) is 

Notwithstanding the fact that, if vp runs to larger values, Og runs to smaller values in 
the infrared, if the fine structure constant is really this large, the accuracy of perturbative 
computations is questionable at best. 

In this paper, we shall formulate a top-down, first principles construction of a class of 
models containing 2+1-dimensional relativistic fermions using the AdS/CFT correspon- 
dence. We believe that these models are as close as one can come with the current state of 
the field to a holographic model of graphene. Moreover, they should have wider applica- 
bility, particularly to topological insulators which can also be strongly interacting systems. 
We will discuss their field theory dual in more detail shortly. For now, we note that they 
and the graphene model (1.1) have some significant differences. For example, the construc- 
tion that we use has 5*0(2, 1) Lorentz invariance and can therefore only describe (1.1) when 
vp = c. In addition, (1.1) has a U{1) gauge symmetry. Our holographic model generalizes 
it to U{N) gauge symmetry, and then analyzes it in the limit A^ — t- oo, a long way from 
A^ = 1. This is not as bad as it sounds. As we will discuss in later sections, if one combines 
the charges of particles gfyM with factors of A^, in the planar limit, A^ itself will make an 
appearance only through the 't Hooft coupling, A = ^^ym-^ which is held fixed in that limit. 
The holographic construction is accurate when this coupling is large and in that sense, the 
large values of planar gauge theory fine structure constant, j;^^ could emulate the large fine 
structure constant Og of graphene. Furthermore, the system which we analyze does have 
gluon-mediated current-current interactions with a -^ propagator very similar in nature 
to the electromagnetic interaction in (1.1). The planar limit of the theory does take this 
into account to all orders in the coupling constant, and the coupling is then taken to be 
large. The main differences are then the addition of higher correlations of the gluon (the 
three and four gluon vertices) and the truncation of Feynman diagrams to those which are 
planar, two effects which compensate each other somewhat. 

Of course there is the question as to whether, if the coupling becomes strong enough, 
a system such as graphene will generate a mass gap for the fermions by spontaneously 
breaking the flavor symmetry (C/(4) in the graphene model (l-l))- There is at present no 
experimental evidence for such a gap, at least in the absence of strong magnetic fields. 
In this paper we will find a strong coupling limit where the fermion spectrum does not 
acquire a gap, but where the strong coupling limit is a nontrivial conformal field theory. 
The different construction, where it does obtain a gap has been examined in references 
[9] and [10]. Our solution at least demonstrates that there is a D-brane construction of a 
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Figure 1. The Feynman diagram with is computed to find the result quoted in eq. (1.6). 

strongly interacting theory where gap generation does not occur. In fact, it turns out to 
be very difficult to find a gapped solution in the class of models that we consider. We do 
find one, but it is only in a system which has extra degrees of freedom. 

A general quantum field theory of the form we want to emulate is described by the 
action 

S = d X ijJkii'j'^da + im)ipk + interactions (1-2) 

Here, and in the rest of this paper, we assume that 2+1-dimensional spacetime has Eu- 
clidean signature. (Dirac gamma matrices are assumed to be hermitian and the Euclidean 
space mass term has a factor of i.) The theory described by (1.2) contains A^ species {k 
is summed from 1 to N) of a 2-component spinor field ipk which are complex. As well as 
gapping the fermion spectrum, the mass term violates parity and time reversal invariance. 
It is absent in the graphene model (1.1) which is parity invariant, but could be present in 
a time reversal violating topological insulator, for example. 
The model (1.2) has a conserved U{1) charge with current 

N 
ja{x) = '^'lpk{x)-fa^k{x) (1.3) 

fc=l 

An interesting probe of the properties of this field theory is the current- current correlation 
function, 

Lorentz invariance implies that it must have the form 



iJaixMO)) = I ^^A,,((7) (1.4) 



Aabiq) = Acs(g) eabcq" + ATiq){q'^Sab - qaqt) (1-5) 

with two functions Acs(9) and A^iq) of the magnitude of q. For A^ two-component Dirac 
spinors with a mass, m, as in (1.2), and at order one loop in perturbation theory, the 
Feynman diagram in figure 1 yields 

N m q 

Acs(q) = 7^ arctan -— (1.6) 

2tt q 2m 






-^ + —(1- ^) arctan -^ 
g2 2q^ ^' 2m 



(1.7) 



These functions have the large and small momentum limits 

N m N 

q»m: Acs{q) = -r- + ■ ■ ■ , Ariq) = -^ + . . . (1.8) 

4 q log 

N m N 

q«m: Acs(g) = ^^ + • • • , Ariq) = -^-^^ + . . . (1.9) 

47r \m\ 12TT\m\ 

respectively. The Fermion mass in (1.2) violates parity and time reversal invariance explic- 
itly and the parity and time reversal violating first term in (1.5) is a result. The leading, 
zero momentum contribution to Acs in (1-9) is independent of the magnitude of the mass 
m and only depends on its sign. This is a manifestation of the parity anomaly [11] and 
leads to an induced Chern Simons term in the effective action of a gauge field which would 
couple to the U{1) current. It also results in a quantum Hall effect, even in the absence of 
magnetic field, with a half-quantized Hall conductivity 

1 e^iV 

Moreover, Acs(9 = 0) has a no-renormalization theorem [12], for a large class of inter- 
actions which preserve the charge gap^, either relativistic or nonrelativistic, corrections 
vanish order by order in perturbation theory beyond one loop. 

On the other hand, it has been demonstrated explicitly that, if the charge gap vanishes, 
i.e. if there are any massless charged particles in the spectrum, Acs(O) can renormalize at 
two and higher loop orders [13]. Then, Acs(O) can differ from that in (1.9). In addition, 
Ar(g) should obtain the cut singularity At ~ l/y/q^ that is expected for massless charged 
particles. The latter is indeed seen at one loop order if one puts m = in (1.2), AT{q) = j^-. 
Both the renormalization of Acs('? = 0) and the nonanalyticity At ~ - are diagnostics 
of the presence of charged massless degrees of freedom. In particular, if a parity and time 
reversal violating conformal field theory existed, conformal invariance and the fact that ja is 
conserved determine the current-current correlator up to constants, where Acs = constant 
and At = constant /g would be exact statements for all momenta. If it were parity or time 
reversal invariant, then Acs(9) would vanish. 

In this paper, we shall find holographic realizations of three scenarios. The first is a 
parity and time reversal invariant conformal field theory which can be regarded as large N 
deformation of the graphene model (1.1) with vp = c. Secondly, we shall find a parity and 
time reversal violating but still gapless theory which flows to a parity and time reversal 
violating conformal field theory at an infrared fixed point. Thirdly, we will find a theory 
with a charge gap whose current-current correlator has functions resembling (1.9) with the 
correctly quantized Acs(O) = -^j consistent with the perturbative no-renormalization the- 
orem. The latter is important. It can be regarded as a nontrivial check of the construction 
- a match between weak and strong coupling of an (albeit coupling constant independent) 



^Here, charge gap refers to the existence of a mass gap for all excitations of the theory which carry 
electric charge. Ordinary 3 + 1-dimensional quantum electrodynamics has no mass gap, since the photon 
is massless, but it does have a charge gap, since the electron has mass. The existence of a charge gap is 
sufficient to make the current-current correlations functions analytic in the region where the 4-momenta 
carried by the currents approach zero. 
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Figure 2. We are studying a defect conformal field theory were fermions are constrained to occupy a 
plane, denoted by the vertical line through the center of the diagram. This plane divides the three 
dimensional space into two different regions which are occupied by four dimensional conformal 
A/" = 4 supersymmetric Yang-Mills theories with different gauge groups. The conformal field theory 
has three tunable parameters, N, N + k and the Yang-Mills coupling constants. The holographic 
description describes the planar limit of this theory where the 't Hooft coupling is tuned to be 
large. The remaining parameter is A: = n^. The anomalous dimensions of operators and the 
current-current correlation functions will turn out to be dependent on this parameter (through / 
in the following sections). 



quantity. All three of our constructions correspond to field theories which at high ener- 
gies approach the parity and time reversal invariant conformal field theory. The parity 
violating conformal field theory we find as the infrared limit of our second solution can be 
regarded as the strong coupling limit of a model of the quantum Hall plateau which uses 
Dirac fermions coupled to a U{1) Chern-Simons gauge theory [14]. Ref. [14] demonstrated 
that such an infrared fixed point exists at weak coupling up to order two loops. Our result 
implies that such a theory can also exist at strong coupling. 

The field theory dual of the system that we shall study is a defect conformal field 
theory consisting of fermions confined to move on a 2-dimensional plane which separates 
3-dimensional space into two regions as depicted in figure 2. The fermions carry a global 
C/(l) charge which couples to the current ja{x) whose correlation functions we discussed in 
eq. (1.4). They also carry the fundamental representation of the gauge group of A/" = 4 su- 
persymmetric Yang- Mills theory and they interact by exchanging gluons and other particles 
in that gauge theory. The interaction-mediating degrees of freedom are allowed to propa- 
gate in the 3-|-l-dimensional space-time whereas the fermions themselves are constrained 
to occupy the 2-|-l-dimensional space-time. The holographic model that we construct has 



Yang- Mills theory with different gauge groups on each side of the defect, as shown in figure 
2. On the D-brane side, this is a result of the fact that, the D7-branes with internal fluxes 
which we shall use are technically D7-branes with nj-) D3-branes dissolved into their world- 
volumes. The D7-brane then forms a boundary between regions with different numbers of 
D3-branes and therefore different amounts of RR flux, thus different ranks of the gauge 
group in the field theory dual. We shall describe the details of our top-down construction 
in the next section. 

The AA = 4 supersymmetric Yang-Mills theory is a four dimensional conformal field 
theory. In particular, this implies that its global color current-current interaction varies 
as J — 2-, there x is the space-time interval, identical to the current-current interaction is 
3+1-dimensional electrodynamics. This feature persists in the holographic construction, 
at strong coupling. One of the mysteries of graphene is the fact that, even though the 
interactions in the model (1.1) are strong, much of the phenomenology is modeled very 
well by free fermions. We shall see that some properties of the conformal field theory, 
like electromagnetic response, that we find in the strong coupling limit bear a strong 
resemblance to what one would find for free fermions. Others, such as the anomalous 
dimension of the mass operator do not, they are tuneable and could differ significantly from 
the their weak coupling limits. As an example of a quantity which is similar, conformal 
invariance fixes the AC conductivity which is a nonzero constant even in the charge neutral 
state of the theory where the electron density vanishes. This behavior is seen in graphene 
and it has been analyzed in the weak coupling regime. Here we see that it is modified very 
little in the strong coupling regime of our model. Experimental estimates of this neutral 
point conductivity vary but are typically of the same order as the results that we find. 

2 Holography 

The holographic construction that we shall pursue has become possible after the devel- 
opment of a series of ideas. First, in order to find a solution of string theory whose low 
energy degrees of freedom are 2-|-l-dimensional fermions, we shall use a system of D-branes 
where the number of Neumann-Dirichlet boundary conditions for the open strings of in- 
terest is #ND = 6. In this configuration, supersymmetry is completely broken. The open 
string spectrum has no tachyon. The only massless states are in the Ramond sector, and 
are therefore fermions [15]. Furthermore, the fermions are chiral in the 2-|-l-dimensional 
sense that they have two-components and their mass term violates the discrete space-time 
symmetries of parity and time reversal invariance. 

This leaves a number of possibilities which are related to each other by T-duality. The 
technically simplest is the D3/D7 configuration outlined below 

lL/ UU UU JU *AJ *AJ dij tJU UU JU 



D2> X X X X 

I?7 X X X 

Db XXX XXX 



(2.1) 
D7 XXX xxxxx 



The N3 D3 and A'^7 D7 branes are extended in 2+1-spacetinie dimensions {x^,x^,x'^) = 
(t,x,y) with 50(2, 1) Lorentz symmetry. The lowest energy states of the 3-7 open strings 
are N^N'/ species of 2+1-dimensional 2-component fermions, similar to those in (1.2).^ As 
a reference system, which will see some use in a later section of this paper, we have also 
included a Z?5-brane oriented so that it forms a #ND = 4 system, i.e. a BPS state, with 
the D3. This is also the effective D5 charge orientation induced on the above L>7-brane 
by an internal flux on the x^-x® plane. Finally, note that the x^ direction is orthogonal to 
both the D3 and D7. Thus, the D3 and D7 can be separated in this direction, introducing 
a bare mass for 3-7 strings. 

The D3/D7 system with this configuration was already used by Rey [9] as a model 
of graphene and by other authors to study the quantum Hall plateau transition [16]. To 
apply holography, the number of Z?3-branes N3 is taken as large and they are replaced by 
the AdS^ X S^ geometry. The L'7-branes are treated as probes and the dynamical problem 
is to find their embedding in AdS^ x S^. The embedding used by Rey is topologically 
equivalent the product space ^^5*4 x S'^ C AdS^ x S^. However, that model is unstable 
to fluctuations of the D7 geometry which violate the Breitenholder-Freedman bound for 
AdSi. The instability has been interpreted as the existence of a phase transition to a 
phase with spontaneously broken chiral symmetry which occurs as the coupling constant 
is increased [9] (see [10] and references therein for recent work on this subject). 

A way to construct a system which is stable at strong coupling was suggested by Myers 
and Wapler [17] who observed that, if one uses the same embedding geometry, but adds 
D7 world-volume U{Nj) gauge fields in a topological instanton configuration on S'^, and 
if the instanton number is large enough, the embedded D7-brane becomes stable, at least 
to small fluctuations of the geometry. An important further refinement of this idea was 
found by Bergman et. al. [18] who observed that, if instead of S"^ C S^, one considers S^ 
as a fibration of two 2-spheres {S'^, S'^) over an interval ip G [0, ^], with metric 

dSgs = dif) + sin il)dsg2 + cos ipds^2 (2-2) 

and the world- volume of D7 wrapping S'^ and S'^, this configuration can be stabilized 
by adding world- volume gauge fields with U{1) Dirac monopole fields on one or both of 
spheres, i.e. turning on a world- volume U{1) flux with 

27ra'Fo = f /f^a + f /f^2 (2.3) 

where (/, /) are constant flux densities and (O2, ^^2) are volume 2- forms on the unit spheres 
(5^, S'^) and where R is the radius of curvature of AdS^. The Dirac quantization condition 
is 

2-Ka'-' N-j ' 27ra''' Nj ^ ' ' 

where ud and ud are integers. This has the technical simplification of using Abelian rather 
than non-Abelian worldsheet gauge fields. As well, it allows one to consider the case of a 
single Z)7-brane, Nj = 1, where the world-volume gauge symmetry is U{1). The result is a 



^The A'' that we used in eqs. (1.2)-(1.9) would be replaced by NzN-j 



model which allowed the authors of Ref. [18] to produce a beautiful proposal for holographic 
duals of quantum Hall states. 

A further word is required about the internal fluxes (2.3). The Z)7-brane action has 
the Wess-Zumino coupling 

^FAFACi (2.5) 



/ 



where C4 is the Ramond-Ramond field in type IIB SUGRA which is sourced by DS-branes. 
Thus a non-zero instanton number density, such as we have in our construction, introduces 
an effective -D3 charge on the Z)7-branes. The fluxes (2.3) induce an integer D3 charge 
AN^ = n£)ho. Since the probe brane carries this charge off to the boundary, this accounts 
for the shift in the rank of the Yang-Mills gauge group across the defect [17] as depicted 
in figure 2. 

In the following, we shall use the construction developed in Ref. [18] with an addi- 
tional ingredient, the requirement that the ultraviolet limit of the theory have the discrete 
symmetries parity (P) and charge conjugation (C). Both of these are putatively violated 
by the construction, C by the background gauge fields Fq and V by the Wess-Zumino term 
in the action for the I?7-brane. Both are restored by augmenting them with additional 
transformations. 

C is easy to fix. It is defined as F — >• —F for all worldsheet gauge field strengths. 
It can be augmented by orientation reversing isometrics of both of {S'^^S'^) so that the 
background field Fq is invariant under the combined transformation. 

The dynamical problem of determining the D7-brane embedding and its world- volume 
gauge fields is obtained from the Dirac-Born-Infeld and Wess-Zumino actions 
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(27ra' 



\2 



Vdet{g + 2TTa'F) - ^ ' iF AFA w^^) 



(2.6) 



where q is the world-volume metric, Nv is the number of D7-branes and T7 = ,„ ,- ,< 
is the D7-brane tension, gg is the closed string coupling constant, to which the radius of 
curvature of AdS^ is related by R^ = AirgsN^ia')'^. The positive overall sign in (2.6) and 
the factor of i in front of the Wess-Zumino term are due to our use of Euclidean signature. 
The Ramond-Ramond 4-form lo^^' is defined by dco^ ' = -ji^^^' where u^^' is the sum of 
volume 5-forms of AdS^ and S^. Generally, a stack of Nf coincident D-branes has a U{N-j) 
gauge field on its unit volume and F in (2.6) is the non-Abelian field strength. Here, we 
are assuming that only the U{1) component of this field strength is excited and the rest of 
the components are set to zero. This assumption is consistent with the equations of motion 
for F. 

We will use the AdS^ x S^ metric 



dS'^ = R^ 



„2 , j„,2 , J,2^ , '^^ , jc<2 

S5 



r\dt^ + dx'^ + dy'^ + dz'^) + ^r- + dS 



(2.7) 



with (iS'es given in (2.2) and 



a;(4) = R\^dt AdxAdyAdz + R'^-^dn2 A (iJ^2 (2. 



with 

d^c^ip) =8sm^ipcos^ilj = l-cosAip (2.9) 

Parity in 2+1-dimensions is the transformation (t,x,y) — )• {t,—x,y) which is an ori- 
entation reversing isometry under which the ^dSs-component of u^^' changes sign. Since 
gauge fields transform covariantly (the background field Fq is invariant), the component of 
F AF on AdS^ also fiips sign. To make the Wess-Zumino term in (2.6) parity invariant, we 
must augment the orientation reversing isometry of AdS^ with a parity reversing isometry 
of 5^, tp ^- ^ — Tp. This transformation interchanges S'^ and S'^. It will be a symmetry 
of the background field Fq only when the fluxes are equal, f = f. For fluctuations of the 
embedding to satisfy the BF bound at large r, it is necessary that /^ > |g [18]. 

We take the ansatz for the L>7-brane which wraps the 2-spheres 5^ and S"^, where 
F = Fq given in (2.3), / = /, it shares coordinates {t,x,y,r) with AdS^ in (2.7) and has 
r-dependent positions z(r) and ^{r), with induced metric 

^ =r\dt^ + dx^ + dy2) + ^(1 + r^z'^ + r^'') + cos^ i^dS]., + sin^ ^jdSl, (2.10) 

where z',ip' = 7/, 7/ • The action and equations of motion are 



S7 = T7 I dr r^ Y^(/2 + 4 sin^ ^) {p + 4 cos^ V) \j\ + r^z'^ + r'^^''^ - fr'^z' 



(2.11) 



(^|:)V ^ L ^ ^Vzr-\f + P.r-4) \ ^ ^ _ 4 cos 2V^ sin 2V^(sin^ 2^ - f) ^^ 12) 



l + (rAV;)2 V V ) '^-'^ V 



,dz if+P^r-yl + {ri^py 
r — = -7= (2-13) 



dr ^ 

V = (/2 + 4cos^V')(/' +4sin4V') - {f + 2|)2 (2.14) 

with T-j = NiTtR^{2tt)'^V2+i where l^+i is the volume of the space spanned by the coor- 
dinates {t,x,y). Here, the fact that z{r) is a cyclic variable in the action (2.11) has been 
used to integrate its equation of motion once and pz is the integration constant. 

In addition to the problem of the embedding of the DT-brane, we shall need to study 
the current-current correlation function in this holographic framework. The U{1) current 
is the field theory dual of the U{1) gauge field whose curvature is F in (2.6). To get the 
2-point function of this current, it is sufficient to study fluctuations of the field strength 
F = Fq + F about Fq in (2.6) to quadratic order, i.e. second order in F. To that order, 
it is consistent to simply insert the solution for the F = Fq geometry into the equation 
of motion for F. For our purposes, it is also sufficient (and consistent with equations of 
motion) to take F to have non-zero components only on the AdS^ space and to only depend 
on the coordinates {r,t,x,y). The action for F at quadratic order is 



Sm = NjT-rilTra'f j d^a ^^det{g + 2'Ka'FQ) g^" g^'F^xFup - i'^F A F A 



J'^ 



In addition to this bulk action, in order to maintain invariance under the RR gauge trans- 
formation which would shift the function c(V') defined in equation (2.8) by a constant, it is 
necessary to add a surface term to this action [18]. The result is equivalent to integrating 
the Wess-Zumino term by parts, and dropping the resulting surface term, so that the action 
becomes 



M) 



(2.15) 
Then, with our Ansatz for the D7 geometry (2.10), and choosing the Ar = gauge, we 
obtain 

Sm = ^^ / d^x j dp (]^{dpAaf + \a^Fl^ + idpc{^)erabcAadbAA (2.16) 

where we have dropped the tildes on F and A^ a, b, c are 2+1-dimensional Euclidean indices, 
Q^^ = (/^ + 4 sin il^){f'^ + 4cos^ ip) and the radial variable has been transformed to 

„^^^ _ /'°°^^ V(/^+4sin^^)(/2+4cos4^) , . 

Note that r = and r = oo are mapped to p = oo and p = 0, respectively. Using the 
fourier transform 

Mx)= f -^e'^-^Aa{q) 
J (27r)2 

the field equation can be written 

[-52 + a2g2 ± 2iqdpc{ij)] A±ip, q) = (2.18) 

and QaAa = 0. We are considering polarization states which obey [igteabc] Ac± = ±qAa±. 
The field equation should be solved with the requirement that the solution is regular at 
the Poincare horizon, p = 00. The on-shell action for the gauge field is then 

S= -^^ j d\\\m[Aa{p,-q)dpAa{p,q)] (2.19) 

The current-current correlator is obtained by taking two derivatives of e~ by the boundary 
value of the gauge field Aa(0, q) and then setting ^4^(0, q) = 0. In the following sections we 
shall quote the results of this procedure for our three solutions of the worldsheet geometry. 
Details of our solution of the Maxwell equation and derivation of current-current correlators 
are given in the appendix. 

3 Parity and time reversal invariant conformal field theory 

The parity invariant solution of (2.12) and (2.13) is 
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which requires that we set the integration constant pz = 0. The D7-brane sits at a constant 
angle and approaches the Poincare horizon at r = 0. The world-volume of the DT-brane 
is AdSi X S"^ X S^, where the 5^'s have radii l/\/2 and the radius of curvature'^ of AdS^ 
is (1 + /■^)/\/l + 2/^. This solution is dual to a conformal field theory which inherits the 
SO {3, 2) symmetry of AdS^. 

The flux / remains as a tuneable parameter. The solution exists independently of the 
number D7-branes N^, as long as A^7 << A''3. The dual field theory has a U{7) global 
symmetry (to apply to graphene, we would set A^7 = 4). 

Note that z = zq when r is infinite and z — )• — oo as r — )■ 0. We interpret this behavior 
as the brane being a fat interface between two regions in 3+1-dimensional space. The 
regions are occupied by theories with different densities, being A/" = 4 Yangs- Mills theories 
with different rank gauge groups. Then, the r-dependence of the solution can be interpreted 
as the change in the apparent position of the brane as the energy resolution is varied. This 
is depicted in figure 3. 

Fluctuations about the solution determine the conformal dimensions of operators. For 
example, a solution of (2.12) which goes to j at large r must have the asymptotic behavior 

V^(r^oo)=f + ^ + ^ + ... (3.2) 

where 



A_|_ can be interpreted as the conformal dimension of a parity violating operator O, with 
correlation function fixed by conformal symmetry, 

,.^/ N.^/ XV const. 

{0{x)0{y)) = — 



|2A+ 
\x — y\ ^ 

At weak coupling, O is identified with tp^p, the fermion mass operator [18], and (O) the 
chiral condensate (V'V')- However, operator mixing makes this interpretation obscure at 
strong coupling. The dimension A+ depends on /. At the bound, /^ = §§, A+ = |. 
When /^ = 2' '-^+ ~ ^' ^^^ classical dimension of tptp. When /^ = 1 the deformation 
becomes marginal, A-|_ = 3. In the interval /^ G (^, 1), A_ > ^^ the unitarity bound for 
an operator in 3 dimensions, and it can also be interpreted as the dimension of an operator 
in a different conformal field theory with ipi being the condensate and '!/'2 the source [19]. 
When ijj = J, a and c{ip) are independent of p, the Maxwell equation (2.18) can 
be solved exactly and the current-current correlator extracted by taking two functional 
derivatives of the on-shell action (2.19) by boundary data Aa{p = 0, q). The result is 

At= ' , Acs = 3.4 

27r^ q 

Acs vanishes due to parity symmetry. The function At((7), whose dependence on q is 
determined by conformal symmetry together with the fact that ja is a conserved current, 



These radii are all in units of R, the radius of curvature of the ambient AdS^ x S 
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r=0 



Figure 3. The conformal field theory of fermions confined to the 2+1-diinensional woldvolume 
of a planar defect, depicted as the vertical curve in the center of the picture, and dividing 3+1- 
dimensional spacetime into two regions. Here, only the coordinate z which is transverse to the 
defect is visible. It is plotted in the horizontal direction. The radial AdS coordinate r is plotted 
in the vertical direction. In the field theory, r corresponds to an energy scale at which one probes 
the system. The region on the left is occupied by A/" = 4 supersymmetric Yang-Mills theory with 
SU{N3 + k) gauge group whereas the region on the right is occupied hy J\f — A supersymmetric 
Yang-Mills theory with SU{N^) gauge group. Here, k — n\, is the number of ZJS-branes that are 
dissolved into the D7-brane. The fact that the interface moves to the left as r decreases from infinity 
is interpreted as the apparent position of the interface changing as one changes the wave-length at 
which it is observed, from zq in the ultraviolet limit (r — )• oo), moving all the way to z = — oo in 
the infrared limited (r — )■ 0). 

now determines the current-current correlator at strong coupling. It is qualitatively similar 
to the weak coupling result for massless fermions, At = ^f^ but with a coefficient that 

depends of /. Recall that / > y fh, and is typically of order one. 

Note that the current-current correlation function is proportional to N-^ which is sup- 
posed to be taken to infinity. This limit can make sense if we remember that we have 
assumed that the quantum of charge in the current is equal to one. However, in the dual 
field theory, the Yang-Mills theory has a coupling constant sym and it would be more 
sensible to normalize the charge by taking into account this coupling constant. Then the 
current-current correlation functions would have the factor N'j, replaced by N^Qy^i = A, the 
't Hooft coupling, which is held fixed in the large N^ limit. It is natural that the current- 
current correlator is proportional to A at weak coupling. The holographic description that 
we have been discussing here is accurate when A is large and our result indicates that the 
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correlation function still is proportional to A in the strong coupling limit. 

The conductivity in this limit is given by (see reference [17] for a detailed discussion) 






±^y/l-^/u;^ 



27r2 

where we have gone back to Lorentzian signature space-time- and q^ = {uj,q). When the 
wave-vector is put to zero, axx approaches a constant. 

4 Parity and time reversal violating solution 

At small r, a solution of (2.12) must have the form 



V' = ^^+cr^ + ... , ^^ = V! + 16iEf^-i (4.1) 

If we search for a solution where ip depends on r with boundary behavior (3.2) and 
(4.1), once one of the three constants {'^pl,^p2,c) is fixed, the other two are determined by 
requiring that the solution is nonsingular. (Solutions at r ~ exist only when pz = 0.) 
Setting ipi to some fixed value corresponds to turning on a source for the operator O 
with conformal dimension A+, and whose expectation value {O) is then proportional to 
the other constant ^2- Since A+ is positive, the operator is relevant and the dual field 
theory is no longer a conformal field theory. However, it will flow to another conformal 
field theory in the infrared, small momentum limit. The equations (2.12) and (2.13) can 
be solved numerically. An example of a solution is depicted in figure 4. 

For the solution in figure 4, where ■0 is not constant, we cannot obtain an exact solution 
of the Maxwell equation (2.18). However, it is easy to solve in the limits where q is large or 
small compared to other dimensional parameters. Here, the only dimensional parameters 
in the problem come from the boundary behavior and we can use, for example, t/j-^ " to 
compare with q. 

When q is large, as is shown in the Appendix, the solution is identical to (3.4). The 
theory has an ultraviolet fixed point which is identical to the parity and time reversal 
invariant conformal field theory that corresponds to the constant ip = \ solution. 

On the other hand, when q is small, we find 



".- ^2 +••• , Acs = -V(/ Vl -P- arccos /) + ... (4.2) 

where corrections are of higher order in q. These functions characterize the electromagnetic 
properties of the infrared fixed point of the field theory dual of the non-constant ^ solution. 
That theory is apparently gapless and parity and time reversal violating. We note that 
Acs(O) differs from the one-loop result (1.9). 

By doubling the degrees of freedom, we could find a parity invariant configuration. 
This could be, for example, an exotic time reversal invariant phase of graphene. In that 
case, two D7 branes behave as we have described above and two D7-branes are their parity 
mirror, with ip{r) replaced by ^ — ij){r). Since the pairs of branes behave differently, the 
global f/(4) symmetry is broken to U{2) x U{2). In this case, this is not spontaneous 
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Figure 4. Numerical solution of (2.12) and (2.13) with / = .8. tp is plotted on the vertical axis, 



i/'(oo) — J and ■0(0) = 5 sin ^ / ~ 0.46. z is plotted on the rear axis and we see that z 
small r and z ^constant at large r. 



1/r for 



breaking, instead it is explicit breaking by turning on the boundary condition (3.2) for two 
of the D7-branes and a similar boundary condition but with ipi,tp2 replaced by —ipi, —^2 
for the other two D7-branes. At weak coupling, this would correspond to turning on the 
parity invariant mass that corresponds to a charge density wave in graphene which was 
discussed in reference [3]. At weak coupling, this would seem to gap the spectrum, but 
for this particular solution, at strong coupling gapless charged excitations seem to survive. 
The parity even part of the current-current correlator At would be as we computed above 
with Nf = 4, however Acs would cancel, and would therefore be zero, consistent with 
parity and time reversal invariance. On the other hand, it would re-appear in a flavor 
current correlation function, and it would give the so-called "valley Hall effect", originally 
described in reference [3]. The value of Acs in (4-2) with Ny set to 4 describes the valley 
Hall effect at strong coupling. 

5 Solution with a charge gap 

It is interesting to ask whether we can find a theory with a charge gap. Charged particles 
are the low energy modes of strings which are suspended between the Poincare horizon and 
the Z)7-brane and the mass of a string state is roughly proportional to its length. In the 
previous two examples, the charged particles were massless, corresponding to the fact that 
the Z)7-brane comes arbitrarily close to the Poincare horizon and the 3-7 strings could be 
arbitrarily short. To get a gapped solution, we need a configuration of Z)7-brane which 
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does not approach the Poincare horizon, something similar to a "Minkowski embedding" 
where the D7 brane pinches off before it reaches the Poincare horizon at r = 0. 

The Z)7-brane can pinch off smoothly when one of the 5^'s collapses, for example, when 
V' — ?• 0. However, this would not be compatible with the existence of the magnetic flux / 
on S"^ unless there is a magnetic source. Here, a magnetic source would be supplied by n^) 
DS-branes where no is the number of units of monopole flux. The flat space configuration 
of the D5-branes is given in the table in (2.1). On AdS^ x S^ , each Z)5 brane wraps S"^ C S^ 
and has nD units of Dirac monopole charge on S'^ . 

However, it can be argued that such a Z)7-brane with a shrinking S"^ attaching to a 
suspended Z?5-brane is not stable, the D7 always has smaller tension than the Db. If the 
Db where connected to the horizon, it would simply pull the D7 brane to the horizon. 
This would seem to rule out the possibility of finding gapped solutions using L>5-branes 
suspended between the Z)7-brane and the Poincare horizon. 

The only alternative is that the 1)5 brane is connected to r — )• c« and that the imbalance 
of tensions pulls the D7-brane back to infinity. We can indeed find numerical solutions of 
the D7 equation which behave in this way. An example is depicted in figure 5. The D7- 
brane begins at {r,ip,z) = (oo, ^, — 0.46). As ip decreases, r decreases until it reaches a 
minimum. Then it starts increasing and returns to the asymptotic region at (r, ^, z) = 
(oo, 0, 0.46). During this interval, the coordinate z increases steadily over a finite range. At 
the latter endpoint, S'^ has collapsed to a point, leaving a source with ud Dirac monopoles 
on the D7 world-volume. We can think of the fiux on the collapsing S'^ as being sourced 
by a D5-brane sitting at r = oo. 

The solution depicted in figure 5 is fundamentally different from the solutions that 
we have found in the previous two Sections since it goes to the boundary at two different 
locations. On the AdS^, these locations are separated in the coordinate z, so they nominally 
correspond to two different dual quantum field theories. 

The operator corresponding to fluctuations of ip, for example, in the asymptotic regime 
where V ~^ f behaves as in equation (4.1), whereas in the regime where V' ~^ it has 
A_ = 1, A-i- = 2, independent of / and identical to the quantities for a D5-brane with flat 
space configuration depicted in (2.1). 

In addition, to solve the Maxwell equation on the world-volume requires two sets of 
asymptotic data. For example, if we require that the world-volume gauge field goes to 
Aa{x) at the V' — ^ f asymptote and Aa{x) at the ^ — )• asymptote, it is straightforward to 
solve the large momentum limit where we obtain two decoupled currents (see the Appendix 
for the details) 

N^Njf + l 2. \ /■~-\ n 

iJaJb) = -^^ [q dab - qaqb) , {jajb) = 



{jaJb) = -^^ \q dab - qaqb) (5.1) 

The j — j correlator reproduces the high energy limit (3.4) of our previous solutions. The 
j —j correlator produces what would be expected for the solution of the L'5-brane geometry 
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pz=l 



Figure 5. Numerical solution of (2.12) and (2.13) with f — .8 and Pz — 1- 4' is plotted on the 
vertical axis. It goes smoothly between {l/r,tp,z) = (0,^,0.46) and (1/r, ?/;,z) = (0,0,-0.46). 

with the constant angle 1^ = 0. The field theory dual to this the D5-brane geometry is well 
known [17]. 

On the other hand, a small momentum expansion is diagonalized by linear combina- 
tions of the currents with 



{j+aj+b) = ,_ f-aMc + 



where 



U-aj- 
Pm = 



vr Pyti 

d^ •y/(/2+4sin^V))(/2+4cos4t/>7 

' mm ^ 

N3N7 r'\.... -----^i_MV 



vr^ 



(i^(l — cos4'(/;) 



Pm J 



(5.2) 
(5.3) 

(5.4) 
(5.5) 



where r^^^ in (5.4) is the minimum value of r that the solution reaches. Corrections in 
(5.2) and (5.3) are of order q^. 

The currents j± are linear combinations of j, j which are normalized so that the charges 
obtained by integrating the time components oi j± are integers, as were the charges of j, j. 
We then observe that the j+ current has the value of Acs(O) that would be expected 
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Figure 6. The field theory which is dual to the gapped solution has a sandwich containing SU{N + 
k) supersymmetric Yang-Mills theory in the interior, sandwiched between two planar defects, one 
corresponding to the D7-branes and the other corresponding to the D5-branes, each with their own 
set of degrees of freedom, and with SU{N) Yang-Mills theory on the exterior. The field theories 
on the planar defects are highly correlated. Here, r is the vertical and z the horizontal directions. 
The width of the sandwich is the range of z which is finite in the ultraviolet, large r regime, but 
shrinks and disappears at a scale which is the minimum of r. 

for a system of N^Nj fermions with a parity violating mass gap. This is identical to the 
non-interacting one-loop result for massive fermions. 

Note that j_ has a superfluid-like pole in its 2-point function, indicating spontaneous 
breaking of a phase symmetry. This breaking is what would be expected from joining of 
the D7 and D5-branes. Figure 6 depicts the geometry of the field theory configuration. 

6 Remarks 



A few remarks about the solutions that we have found are in order. First, if we begin with 
the Z)7-brane embedding equation (2.12), which we recopy here for the reader's convenience, 



(^|:)V 



+ 3 + 



Ap^.r ^{f'^ +Pzr~ 



^ 4cos2ipsm2'tp{sm'^ 2ip — p) 



l + irf,^? V y y .r. y 

V = (/2 + 4cos^V')(/' + 4sin^V') - {f + ff)^ 
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We recall that, to find a solution that extends to the Poincare horizon at r = 0, we must 
set pz = 0. We can re- write this equation as 



r— - in 
dr 



(/^ + 4 cos^ V^) (/^ + 4 sin^ ip) - f^ 



d ^' 



6 ( r-V^ ) (6.1) 



The right-hand-side is positive and we can conclude that the left-hand-side is a monotoni- 
cally increasing function of r. 

If we assume that the logarithmic derivative of "0 vanishes at both limits, r — )■ oo and 
r — )• 0, integrating this equation yields the sum rule 

if + 4cos4 V(oo))(/2 + 4sinV(oo)) -f_ ^^^ / f^ dr (^±y\ (g 2) 



(/2 + 4cos4V'(0))(/2+4sin4V(0))-/^ \ Jo r \ dr 

This indicates that any solution of the equation of motion will necessarily have smaller 
V = (/2 + 4sin^'0)(/2 + 4cos^^) — /^ at r = 0, that is in the infrared limit, than at 
r = 00, the ultraviolet limit. If we interpret the evolution from large r to small r as a 
renormalization group flow, we have found a quantity which definitely decreases. What is 
more, it is directly related to the Routhian'* 



^^ ^ N3NjV2+i r ^^ ry(/2 + 4sin4 i,){p + ^cos^i;) - pJl + (r-^A (6-3) 

whose variation gives the equation of motion. It is the density in this integral which must 
decrease, an analog of an H-theorem which has recently been discussed for nonsupersym- 
metric 2+1-dimensional field theories [20]. 

If we evaluate the left-hand-side of (6.4) with the solution that we have found, it reads 

The left-hand-side is greater than one for all values of /^ < 1 and it equals one when 
/2 = 1. When / — )• 1, the operator that we have perturbed the conformal field theory by 
when we switched on ?/^i approaches a marginal operator. It should be possible to develop 
a perturbative approach where this flow can be analyzed explicitly. 

As another observation, in the same solution, the scale symmetry of the embedding 
equation for the L)7-brane is broken only by one of the parameters in the boundary condi- 
tion, say ipi. Then the other parameters in the boundary conditions are determined by ipi 
and the relationship is governed by scale symmetry, 

V'2 = 5(/)V'f ■'/^- 
It would be interesting to compute the function g{f) to see if it has any special behavior. 



^This is obtained by the Legendre transform R{ip{r),'ip' {r),pz) — L{ip{r),ip'{r), z'{r)) — z'{r)pz and 
setting pz = to focus on embeddings entering the Poincare horizon. 
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Finally, the gapped solution that we have found in section 5 seems to be the unique 
example of a simple solution with a mass gap. We have outlined why it is unlikely that a 
solution could exist with the Z)5-brane reaching the Poincare horizon, rather than being 
located at the large r boundary. We believe that the argument is robust. We cannot rule 
out more complicated complexes of branes with multiple intersections. The solution that 
we did find had an additional asymptotic region which essentially doubles the degrees of 
freedom. It also has a spontaneously broken U{1) symmetry which we find explicitly when 
we study the current-current correlation functions. This breaking could proceed with a 
■0^-condensate which breaks the U{1) x f/(l) gauge symmetry to a diagonal U{1). This 
would be sufficient to gap the spectrum of fermions which live on both the D7 and D5 
branes. 
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A Solution of the world-volume gauge field 

The equation of motion for the worldsheet gauge field A± is given in equation (2.18) as 

[-92 + q^a\p) ± tqdp{2c{p))] A± = 

In this Appendix we will outline our solution of this equation for the three different world- 
sheet geometries that we have discussed, beginning with the simplest one with constant 
angle ■(/'• 

Constant a and c: In this case, we have 

[-a^ + q^a^] Aa = 
where a = f'^ + 1. The general solution is 

A = cie'^'^P + 026"''"'' 

The solution must converge at p — )■ 00 (this is r — )■ 0, the Poincare horizon). To satisfy 
this condition, we set ci = 0. Then the solution is 

Aa{p,q)=Aa{q)e-''"f 
Now, we must plug this solution into the on-shell action (2.19) to get 

S = ^ ^ / S'qAa^-q) -{q^Sab - qaqb)Ah{q) 
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Taking two functional derivatives by Aa{q) yields 

Acs = 

iVaiV? f + 1 

which is the result quoted in equation (3.4). 

Non-constant ip and large \q\: Consider the Maxwell equation at large q where now a 
depends on p, 

[-92 + gV(p)]A± = 

and consider the WKB Ansatz 

A = Ao{q)e'^ 

so that the differential equation becomes 

~W"-{W'f + q^a\p)=0 

For large q, W ^ q and 

fp 
W = — dpqa(p) 
Jo 

We have chosen the minus sign so that the solution converges at /o — )• oo. The on-shell 

action is 

S = ^ I d'qW'Ao{q)Ao{-q) = ^ J d^qA,{-q)^-^^j^{q'6at - qaqb)AM 

This gives At, Acs identical to the constant solution that we analyzed in the above sub- 
section (as we expected). 

Non-constant -0 and small q: We can re-scale p in the Maxwell equation by a factor 
of g: p — 7- p/q to obtain 

[-dl + a\plq) + id,Clc{plq))\ A+ = 

and a similar equation for A- which is the complex conjugation. 

From the numerical solution in figure 4, we see that the interpolation of ?/^(r) between 
its two asymptotic regions occurs in a relatively small interval. In the small q limit, this 
interval is squeezed into the small p region. To deal with this, we approximate a and c by 
step functions at a given value of po- Then we will take the limit as po goes to zero. 

In the p < pq regime, a = (/^ + 1)^. In the p > po region, q^ = 4/^ and we approximate 
a^ by the step function 

a2 = (;2 ^ -^)2^(^^ _ ^) ^ i2ffeip - po) 
We will approximate the term with a derivative of c{tp) by a delta function. Since 

i arcsin / i i 



c(oo) - c(0) = /_ diid^c = - arccos / - -/\/l - P 
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we write 



f-fVl^]sip-po) 



9p(2c('0)) = arccos 
The solution of the equation is then 

A+ = [ylie^-^'+i)^ + A2e-^^'+M e{po - p) + Be-^^Pd{p - po) 

where we have imposed convergence at large p. Continuity of the wavefunction at p = po 
implies 

and the discontinuity of the derivative due to the delta-function potential requires 
if + 1) Uie^^'+i)"" - ylae-^^'+i)""] - 2fBe-'fP° 



arccos 



/ _ /yrrp Aie(^ +i)^« + ^2e-(^ +1)''° 



From these equations, we can obtain 



Ao 



2/ . A/1 -/2- arccos/ 



P + 1 



P + l 



{A1+A2) 



Then, the on-shell action is 



S 



47r2 p-)-o 



mNj P + l 
limJA9,A) = - ^J^ , (Ai + ^2)(Ai - A2) 



27r2 



= / + i(/yr^72- arccos /)J (^1+^2)' 
The current-current correlator components are 

N^Nj 2/ 



A, 



cs 



At 

N■iN^ 
1^ 



27r2 



(/\/l - P - arccos /) 



as quoted in eq. (4.2). 



Solution with a charge gap at large q: We use a variable s = p for the interval 
/9 G [0, pm\ and s = 2pm — p for p £ [pm, 0]. The differential equation is 

[-92 + q^a\s) ± iqds{2c{^{s)))] A± = 

There are two boundaries, one at s = 0, the other at s = 2pm- Let us assume that 
A{0,q) = A{q) and A{2pm,q) = A{q). 

We first study the large q'^ regime. The solution is 



A{s) 



-Ae-'^+Ae-^ A-Ae 

e^ + 



-w 



-w 



l-e 



-2W 



l-e 



-2W 
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where W[s] = Jq ds'a{s')q and W = W[2pm]- Since W ^ q we should take the hmit where 
W is large. 

d,A{0) = -{f + 1) coth WA ~ -{f + l)A 

dsA{2pm) = fyJP + AcoihWA ~ f^/p + A A 

Then the on-shell action is 



S 



NsNj 
~2^ 



fVP+^ 7 . 2 e w~ ^ 
Aa{q dab - qaqhjAb + 

Iq 



(f + 1) . .„2 



2q 



Aa{q 6ab - qaqb)Ab 



which reveals the result quoted in eq. (5.1). 



Solution with a charge gap at small q: First, we observe that, when g — )• the 
Maxwell equation is solved by 



A{s) = A{q){2p^ - s)/2pmax + s/2pmA{q) 
and the on-shell action is proportional to 



1 



Ad, A 



s=2p„ 



1 



4p„ 



■\Aiq) - Aiq)(^ 



This is a Higgs-like term, it implies that the difference of currents j — j gets a pole in its 
two point correlator. 

Now, to expand the solution in q, it is convenient to convert the differential equation 
to an integral equation, 

A = A{q) (l - -^) + -^A{q) + r ds' r ds" [q^a\s") ± iqds"{2c{s"))\ A{s") 

\ ^PraJ ^Pm JO JQ 

j'2pm fs' 

--— / ds' / ds" [q'^a'^is")±iqds"{2c{s"))] A{s") 

^Pm Jo Jo 

and then we can solve iteratively. We shall only examine the solution to leading order in 
q. Iteration yields 



A{s) = ^ ( 1 - -^ ) + -^A ±iq ds' / ds"d,.{2c{s")) 

2pm/ 2pm Jo Jo 



A 1 



s" \ s" 



+ A 

■^Pm / ^Pm 



2pm rs 



Then 



=fiq- / ds' / ds"ds"{2c{s")) 

2pm Jo Jo 



2pn 



A 1 



" \ J/ 



■^Pin / -^Pn 



-A 



1 _ / "f"- a' 

d^AiO) = ,— {A -A)^iq ds'{l - ,—)dA2c{s')) 

^Pm Jo ^Pm 



A I 



+ 



+ z^A 



2pm J 2p, 



+ ... 



dsA{2p„ 



1 



2pm. g' 



-{A-A)±iq ds',—dAHs')) 

'^Pm Jo ^Pm 



^n.^U^i 



^Pra J ^Pn 



+ ... 
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The on-shell action is proportional to 

1 



AdsA{2pm)-A{Q,)dsA{id) 



2pri 



{A-Af±iq I ds' ds'{2c{s')) 

^Pm Jo 



A 1 



2pn 



The matrix of current-current correlators is 



< jj > < n > 

< jj > < jj > 



2^2 



Om ^ ^ Pm 

--j-±iqk2 -j-±iqks 



— ±iqki -—±iqk2 

Pm Pi" 



where 



2/ ds' dAHs'))(i-,^ 

/o V 2pm 



2p„ 



„/ 



k2 = 2 I ds' dA2c{s'))[l ^ .^ 

Jo \ -^Pm/ ^Pn 

r2pm 

k3 = 2 ds' ds'{2c{s')) . ^ 

Jo \ ^Pm 

To linear order in q, the eigenvalues are 

^ . . ^ , . N^jNj {ki + 2k2 + ks) 

< J+J+ >= ±^q^:^2 2 + • • • 

/ . . ^ 2 NsNr (ki - 2k2 + ks) 

< 3-3- >= ± ^1^r2 n + ■ 

Pm 27r^ 2 

N3N7 

< J+J+ >= T«g— ; h . . . 



or 



< j-j- >-- 



Air 

± ^Q n 9 / ds ds'{2c{s )) 1 - 2- — 
~ /o V 2p„ 



+ 



2p„ 



-A 



+. 



+ ... 



2vr2 p,„ 27r2 

which are the results quoted in eq. (5.3) and (5.2). Furthermore, in this approximation, 
j± = j ±j + 0{q) and the elementary charge quanta for j+ and j- are identical to those 
for j and j. By convention, we have taken these charges to be integers. 
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